Superconducting transitions from the pseudogap state: rf-wave symmetry, lattice, 

and low-dimensional effects 
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We investigate the behavior of the superconducting transition temperature within a previously developed BCS- 
Bose Einstein crossover picture. This picture, based on a decoupling scheme of Kadanoff and Martin, further 
extended by Patton, can be used to derive a simple form for the superconducting transition temperature in the 
presence of a pseudogap. We extend previous work which addressed the case of s-wave pairing in jellium, to 
explore the solutions for Tc as a function of variable coupling in more physically relevant situations. We thereby 
ascertain the effects of reduced dimensionality, periodic lattices and a d-wave pairing interaction. Implications 
for the cuprate superconductors are discussed. 



PACS numbers: 74.20.-z, 74.20.Fg, 74.25.Dw, 74.25.Jb 



I. INTRODUCTION 

The concept of a smooth evolution from a BCS descrip- 
tion of superconductivity to thatj^f Bose Einsteio-condensa- 
tion (BEC) dates back to EaglesEl and to Leggett.El The latter 
addressed this problem at zero temperature in the context of 
p-wave pairing in He^. Nozieres and Schmitt-Rinku (NSR) 
extended Leggett's formalism to calculations of Tc and, for 
the case of a jellium gas, found a continuous variation from 
the BCS exponential dependence (on coupling constant g) 
to the Bose-Einstein asvrnptote, at large g. UemuraQ and, 
independently, RanderiaQ and MicnasQ and their respective 
co-workers applied this BCS Bose-Einstein crossover picture 
to the high-temperature superconductors, which, because of 
their short coherence length, were claimed to correspond to 
intermediate values of the coupling. It was, subsequently, 
argued by these and other groups that the pseudogap (nor- 
mal) state of the cuprates was naturally associated with this 
intermjediate coupling regime. Since then, a large number of 
paperstrEl have been written on the crossover problem and the 
related pseudogap state. 

The application of these crossover theories to the cuprates 
is made complex by a number of important factors which 
involve quasi-two dimensionality, lattice periodicity, and, fi- 
nally, the introduction of c?-wave symmetry in the pairing in- 
teraction. It is the goal of the present paper to discuss these 
three effects in the context of a many-body-theoretic approach 
to the crosspijer problem, based on earlier work by Kadanoff 
and Martinjl3 and extended by Patton.U A major advantage 
of this scheme is that it reduces to BCS thepB|!4n the limit of 
weak coupling. Our principal contributionsll3ll^ to the larger 
body of work on the BCS-Bose Einstein crossover have been 
based on s-wave pairing in three-dimensional (3d) jellium. 
In the context of this diagrammatic approach we have es- 
tablished that (i) there is a breakdown (for T < T*) of the 
Fermi liquid at intermediate g, which roughly coincides with 
the onset of long lived (i.e., resonant) pairs and (ii) that this 
breakdown has characteristic pseudogap features, such as a 
depressed density of states at Fermi energy Ep, as well as a 
two-peaked spectral function. The pseudogap amplitude Apg 
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can, moreover, be quantified, (iii) We have examined the su- 
perconducting instability associated with this pseudogap state 
and determined the three self-consistent, coupled equations 
which must be satisfied for the inter-related quantities Tc, ^, 
and Apg(Tc). (iv) In the process we have presented a quantita- 
tive phase diagram for Tc and T* as a function of the coupling 
9- 

Various aspects of low dimensionality, lattice effects, and 
non-s-wave pairing interactions have been addressed in the lit- 
erature withiiLa crossover scenario. Schmitt-Rink, Varma, and 
Ruckensteiro applied the NSR approach to two-dimensional 
(2d) systems and found a breakdown of the Fermi liquid even 
for arbitrarily weak coupling g. This was manifested as a 
negative chemical potential ji, which occurred in conjunction 
with Tc = 0. Serenella suggested that this breakdown wp. 
an artifact of the NSRrScheme, which is not conserving.E3 
Yamada and co-workersllZl introduced a diagrammatic "mode- 
mode coupling" scheme, following similar work on related 
magnetic problems, and they found that /i was properly pos- 
itive at weak coupling, while Tc remained zero, as expected. 
However, they were unable to find a continuous crossover be- 
tween the limits of very strong and weak interactions. 

Lattice effects were discussed in the original Nozieres pnd 
Schmitt-Rink paper (and later by Belkhir and Randerialla). 
These authors noted that the Tc calculations of the jellium 
case, could not be readily extended to the lattice, at least in 
the strong coupling regime. This difficulty came from a vari- 
ety of issues, among which, it was claimed, was the neglect 
of interactions between pairs of composite fermions. In ad- 
dition there is a reduction in the effective pair hopping matrix 
element or kinetic energy. These combined effects conspire so 
that Tc is expected to vanish at arbitrarily large g, in contrast to 
the transition in the Bose-Einstein ideal gas. Demonstrations 
of this effect, in the negativej/ Hubbard model, came later, 
through Monte Carlo studies,Ej as well as from additional an- 
alytic work based on the coherent potential approximation.^ 

The role of non-s-wave pairing symmetry was initially iiu. 
vestigated by Leggett and, subsequently, by otlier groupsEll 
which addressed jellium models. It was noteda that in the 
strong coupling limit the effects of an anisotropic order pa- 
rameter are absent in the excitation spectrum, so that there 
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can be no gap nodes in the bosonic regime. t/-wave pairing 
on a lattice is not yet amenable to Monte Carlo approaches 
because of the fermion sign problem. There have been, nev- 
ertheless, some numerical solutions which address electronic 
spectral functions, based on the "fluctuation epaphange approx- 
imation" (FLEX) for the lattice d-wave case.tj Because these 
were applied to a strictly two dimensional model there was no 
discussion of the behavior of Tc in the more general context 
of the crossover problem. 

The results of the present paper, in which Tc is computed as 
a function of g, can be put into the context of this background 
literature. Here we arrive at a crossover scheme which in the 
strict 2d limit yields Tc — and in which /i smoothly interpo- 
lates from Ep to large negative values in the jellium case, as 
the coupling varies from weak to strong. Lattice effects yield 
a vanishing strong coupling limit for Tc, associated with a re- 
duction in the effective kinetic energy of the bosons or pairs. 
The c?-wave case on a lattice is found to be different from the 
s-wave lattice case in one significant respect: we find that su- 
perconductivity disappears at relatively smaller values of g, 
as a result of d-wave symmetry; the pair size cannot be less 
than a lattice spacing so that the pairs interact more strongly. 
As a consequence their mobility is suppressed. Consequently, 
the bpspnic regime is essentially never reached for the d-wave 
case.E2l The implications of this observation and other aspects 
of these calculations for the underdoped cuprates are briefly 
addressed at the end of this paper 

II. THEORETICAL FRAMEWORK 

We consider a generic system of fermions characterized by 
an effective, short range pairing interaction with Hamiltonian 

kCT 

+ ^ ^k,k'4+q/2TC-k+q/2iC-k'+q/2iCk'+q/2T' (1) 
kk'q 

where cj^^ creates a particle in the momentum state k with 
spin a, and ek is the energy dispersion measured from the 
chemical potential /i (we take h = ks — 1). For simplicity, 
we assume a separable pairing interaction Vk.k' — ffVk'/'ks 
where g — —\g\ is, the coupling strength; the momentum 
dependence of the function (^k, which reflects the pairing 
anisotropy, will be specified below. 

In order to establish notation and to make clear our approxi- 
mations, in what follows we provide a brief description of our 
formalism. This approach is based on the correlation func- 
tion (or Greens function equation of motion) formulation of 
the "pairing-approximation" originally discusjsed by Kadanoff 
and MartinEJ and later extended by Patton.EJ The system can 
be characterized by the one- and two-particle Green's fLmc- 
tions which, in real space, obey the following equations:t3 

G(l - 1') = Go(l - 1') (2a) 

+ /dTd2Go(l-T)1^(T-2)G2(T2;l'2+) , 



G2(12;1'2') = - Jdld2G{l - l)Go(2 - 2) (2b) 

xV(l-2) G2(T2;1'2') , 

where the two-particle correlation function G2 is given by 

G2(12;1'2') = G2(12;l'2')-G(l~l')Go(2-2') 

(2c) 

and, for brevity, we have used a four vector notation 1 = 
(r,T), etc. While Eq. ( ^ ) is exact, Eqs. ( ^b| ) and (^ are 
approximate; these eq-uations were originally proposed by 
Kadanoff and MartinEiLJ [see their Eqs. (2.6) and, most im- 
portantly, (2.29)] as a simple decoupling scheme for the three- 
particle Green's function. 

This approach is primarily motivated by the observation 
that it yields the results of BCS theory, in the weak coupling 
limit. It is convenient to express the correlation function G2 
in terms of a T-matrix (or pair propagator) via the definition 

G2(12;l'2')= Jdld2dl' d2' G{1 -T)Go{2-2) (3) 

xt(T2;T'2')G(T'-l')Go(2'-2') . 

Taking the Fourier transform of Eqs. ^ and (|^), and noticing 
that for our separable interaction the T-matrix can be written 
as t{K,K';Q) = t{Q)ip]f^ipi^i, after sorae. straightforward al- 
gebra we obtain the following equationsO for the self-energy 

= 5^t(Q)Go(g~X)^Lq/2: (4a) 


and the T-matrix 

9 = [l + 9XmtiQ) , (4b) 

where 

xiQ) - 5]G(i^)Go(Q-if)^Lq/2 (4c) 

K 

is the pair susceptibility. Here, Go{K) = l/{iuj — ek) is the 
bare propagator, and K = (k, icj), J2k = ^J2k.i<^' ^tc, 
with n/ui denoting the even/odd Matsubara frequencies. To- 
gether with the particle number equation 

n = 2^G(if) , (4d) 

K 

Eqs. (0) form a complete set which, for a given g and T, need 
to be solved self-consistently for 'S{K), t{Q) and the chemi- 
cal potential p.. Equations (Q) represent the basis of the present 
theory and can be regarded as a generalized BCS theory in 
which pair correlations are explicitly taken into account in the 
normal state (T > Tc) in a self-consistent fashion. 

This approach has implications for the superconducting 
state as well. The presence of (noncritical) pair fluctuations 
in the supeconducting state (T < Tc) leads to progressively 
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stronger deviations from BCS theory as the coupHng is in- 
creased. However, this standard (BCS) theory is embedded in 
Eqs. (0) when the T-matrix is given by 







S{Q) 



for r>T,, 
for T<T,, 



(5) 



where Asc is the superconducting order parameter. Insert- 
ing Eq. into Eqs. one obtains (i) the usual BCS self- 
energy J:bcs{K) = \A^c\'^ipl/{iuj + e_k) and (ii) gap equa- 
tion = 1+ .9 x(0) = 1 + g Sjf |Asc|Vk/(^' + El), 
where i?k = -^/e^ + | Ascpi^k '■^^ energy dispersion of 
the quasiparticles. The delta function in t^dQ) leads to 
the usual (Gor'kov) factorization of the correlation function 
C2{K,K') = F{K)F{K'), where the anomalous Green's 
function (if) = Asc(pkGoi-K)G{K). 

Under more general circumstances, when pair fluctuations 
cannot be neglected, the self-consistent T-matrix in Eqs. (Q) 
can be written as 



t{Q) 
ipg(Q) 



tsc(Q) + tpg(Q) 
9 

l + 9XiQ) ' 



(6a) 
(6b) 



where tpJQ) is the "regular" (or pseudogap) contribution 
[cf. Eq. (ftb|)], which should be associated with (noncritical) 
pair fluctuations which persist both above and below Tc- At 
the high temperatures of the normal state, tpg{Q) is finite at all 
Q. As the temperature is lowered tpg(Q) develops a resonant 
structure corresponding to metastable or long lived pairs. Pre- 
cisely at Tc, this quantity becomes divergent for Q = 0, in ac- 
cord with the Thouless criterion for a superconducting pairing 
instability. Once the temperature is less than Tc, a nonzero su- 
perconducting order parameter Asc is established which obeys 
the (gap) equation 1 + 5x(0; T) = 0. It is important to stress 
that the same critical temperature is obtained either when ap- 
proached from the normal state (using the Thouless criterion) 



C(0)=.9"'+X(0;T,)=0 



(7) 



or when approached within the superconducting state by set- 
ting Asc — within the gap equation. 

The calculation of Tc can be substantially simplified by not- 
ing that atnthis temperature the self-energy is well approxi- 
mated bycJ 



(Q) 



(8a) 



This approximation is highly nontrivial and establishing its 
validity requires detailed numerical calculations which are 
presented elsewhere.Ej It should be stressed that Eq. ( p^ can- 
not be written down on analytical grounds-alone. It is a con- 
sequence of a numerical iterative solutionllj of Eqs. (Q) and 
obtains when the T-matrix contains a divergence. In the pseu- 
dogap phase Eq. ( p^ ) is no longer a valid approximation to 
Eq. (p^. Indeed, one may arrive at some apparent inconsis- 
tencies if this approximation is used above Tc- It follows from 
Eq. (pa) that at Tc, the self-energy has the same symmetry as 



in the superconducting state. However, it can be seen from the 
exact Eq. (^, that the pairing symmetry factor (pj. cannot, in 
general, be extracted from inside the integral. When the in- 
tegration is properly performed, the anisotropy at higher tem- 
perature will generally be different from that at Tc- This im=. 
portant consequence, which may be relevant to experiment£3 
would not follow if Eq. (|^ were incorrectly extended beyond 
its regime of validity. 

The self-energy can be rewritten as 



S(if) 



(8b) 



where the pseudogap parameter is defined as 

/CXD ir-\ 
—b{n)lmt^,n (9a) 

in terms of the Bose function b{^l). Here, and in what fol- 
lows, we use the notation Apg to represent the amplitude of the 
pseudogap at Tc, while the momentum dependence is given by 
Ak = Apgi^k- The simple BCS-like form of the self-energy 
( |8b| ) allows us to express Eqs. (^ and ( |4d| ) as 



and 



2E 



Ck 

Sk 



(9b) 



(9c) 



where v\ = i(l - ek/Ei^),u\ = i(l + CkZ-Ek), and 
f{E) is the Fermi function. The complete set of Eqs. (^), 
which were previously established in a slightly different form 
in Ref. ^ must be solved self consistently in order to obtain 
Tc, /i, and Apg as a function of g and n. 

The imaginary component of the T-matrix which appears 
in the pseudogap equation ( pa]), can be directly calculated by 
inserting Eq. ( pb[ ) into Eqs,3^3) and (pc[). However, exten- 
sive numerical calculationsll3 EZTshow that as Tc is approached 
from above, one can simplify this procedure considerably and, 
at the same time, obtain considerable physical insight. For the 
purposes of calculating quantities such as Apg (which involve 
integrals over the T-matrix), it suffices to approximate the T- 
matrix near Tc by its values at small fl and q. We derive the 
appropriate form by noting that the inverse of the analytically 
continued function can be written as 



Ret 
Im< 



q,0 

-1 
q,f2 



a'o {ft - n^) 



where fiq can be naturally interpreted as an energy disper- 
sion of pairs of fermions. Here the parameters a'g and Oq are 
essentially constant iiJ4he relevant range of momentum and 
frequency. MoreoverJ13 independent of the values of g, suf- 
ficiently close to Tc, we find that the ratio e = a'^/a'^^ ^ 1. 
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Hence the imaginary part of the T-matrix can be approximated 

as 



Im t, 



— lim 



Imt 



q,0 



Ret 



q,0 



Imt 



q,0 



lim 



Finally, it is useful to rewrite the approximated T-matrix at 
Tc (for small momenta and frequencies) in a more compact 
form as 



t. 



q,n 



^0 



(10) 



This approximated T-matrix takes the natural form of a pair 
"Green's function" or propagator, with characteristic disper- 
sion riq. Using the inversion symmetry of the Hamiltonian 
we deduce that flq varies quadratically at Tc with the wave 
vector, as 



(i) 3d jellium, s-wave symmetry. We assume a parabolic 
dispersion relation, Ck — k^/2m — /i, with if^ — {1 + 
fc^/fcQ)^^/^. The parameter fco is the inverse range of the 
interaction and represents a soft cutoff in momentum space 
for the interaction. As will be clear later, fco > fc_F is as- 
sumed in general in order to access the strong coupling limit. 
It is convenient to introduce a dimensionless scale g/gc for 
the coupling constant. Here, following Ref. |], we choose 
gc = —iTr/mko, which corresponds to the critical value of 
the coupling above which bound pairs are formed in vacuum. 

(ii) 2d jellium, s-wave symmetry. For 2d jellium we choose 
the same et and (pk as for case (i). We find that Tc = 0, in 
agreement with the Mermin- Wagner theorem. To understand 
this result, note that the assumption that Tc is finite leads to a 
contradiction, associated with an unphysical divergence in the 
pseudogap amplitude. This unphysical result derives from an 
infrared, logarithmic divergence in the phase space integral on 
the right hand side of Eq. (^. This divergence can be made 
obvious by rewriting this equation using the low frequency, 
long wavelength expansion of the T-matrix so that 



(14) 



2Af* 



2M 



- + 



2M-1 



for 3d, 



for quasi-2d, 



(11) 



where M* is the effective mass of the pairs. By quasi-2d we 
mean a highly anisotropic 3d system with 7\/J^/il/|* ^ 1. The 
effective pair mass is determined by an expansion of the pair 
dispersion, given by 



1 - f{Ek) - /(ek-q) 



/(^k-q)^ 



Ek — Ek-q 



Ck-q 



(12) 



where 



2A2 ^ 



[1 - 2/(ek)] 



Ck 

Ek 



[1 - 2/(£;k)] 



(13) 



We have verified numerically that the above leading order (in 
q) contributions in Eq. ( pi] ) dominate in our subsequent cal- 
culations, so that higher order terms in the expansion can be 
dropped. The relationship between the effec tive m ass of the 
pairs and Tc will be explored in detail in Sec. [II A. 

The above analysis will be applied to isotropic and 
anisotropic jellium with s-wave pairing, as well as to discrete 
lattices. In the latter case we consider both s- and d-wave 
symmetry of the pairing interaction I4;,k' • The distinction be- 
tween these various situations enters via the dispersion rela- 
tion ek and the symmetry factor ip-k^, which will be charac- 
terized below, according to the details of the physical system. 
For definiteness, in our quasi-2d calculations it is assumed that 
the pairing interaction depends only on the in-plane momenta. 



Pairing fluctuations, thus, disorder the system for any finite 
temperature. Even in 2d, for which Tc = 0, we obtain a finite 



pseudogap, as will be seen in Sec. Ill B It should be noted that 
this result is general and remains valid for both s- and d-wave 
pairing on discrete lattices, as well. This is a consequence 
of the fact that the lattice energy dispersion is quadratic at 
sufficiently small wavevectors, so that the same arguments as 
above can be applied. 

(iii) Quasi-2d jellium, s-wave symmetry. Here we use (pk, 
as in the previous two cases, and adopt an anisotropic energy 
dispersion 



Ck = 



11 



2mii 



2m± 



(15) 



where k±_ is restricted to a finite interval (|fc_L| < 7r),§l while 
k|| is unconstrained. 

By tuning the value of the anisotropy ratio m±/m\^ from 
one to infinity, this model can be applied to study effects as- 
sociated with continuously varying dimensionality from 3d to 
2d.Ej For convenience, we use the parameter gc derived for 3d 
jellium, as a scale factor for the coupling strength, and call it 
go to avoid confusion. 

(iv) Quasi-2d lattice, s- and d-wave symmetry. In the pres- 
ence of a lattice we will adopt a simple tight-binding model 
with dispersion 

Ek = 2f||(2 — COSfca; — COS ky) +2t±{l— cosfcj_) — fl , 

(16) 

where t|| {t±) is the hopping integral for the in-plane (out-of- 
plane) motion. Here we consider both isotropic s-wave pair- 
ing symmetry with (/?k = 1, as in the negative U Hubbard 
model, as well as d-wave effects with 
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ifk = COS kx — cos ky . (17) 

It should be noted that in the lattice case, because the momen- 
tum integration is restricted to the first Brillouin zone, it is not 
necessary to introduce a cutoff for the interaction in momen- 
tum space. 

m. NUMERICAL RESULTS 

Equations (^, together with the various models for (/3k and 
eic, were solved numerically for Apg, /i, and Tc- The numeri- 
cally obtained solutions satisfy the appropriate equations with 
an accuracy higher than 10"^. The momentum summations 
were calculated through numerical integration over the whole 
k space for the jellium case, and over the entire Brillouin zone 
for the lattice. However, to facilitate our calculations in the 
case of the quasi-2d lattice with a d-wave pairing interaction, 
the momentum integral along the out-of-plane direction was 
generally replaced by summation on a lattice with N± — 16 
sites. For completeness we compared solutions obtained with 
and without the low frequency, long wavelength expansions 
of the T-matrix discussed above, and found extremely good 
agreement between the two different approaches. In general, 
we chose the ratios tox/to|I = 100oit±/t^^ = 0.01, although 
higher values of the anisotropy were used for illustrative pur- 
poses in some cases. 

A. Overview: Tc and effective mass of the pairs 

In was pointed out in Ref. ^ in the context of the attractive 
Hubbard calculations, that the appropriate description of the 
strong coupling limit corresponds to interacting bosons on a 
lattice with effective hopping integral t' « —2t^/U. It, there- 
fore, will necessarily vanish in the strong coupling limit, as 
U ^ oo. In addition to this hopping, there is an effective 
boson-boson repulsion which also varies as V' ~ —2t^/U. 

This description of a boson Hamiltonian can be related to 
the present calculations through Eqs. (|lO|) - (|2|) which rep- 
resent the Green's function for such a Hamiltonian and its 
parametrization via the pair mass M*. By solving Eqs. ^ 
self-consistently and identifying M* from the effective pair 
propagator (or T-matrix), our M* necessarily incorporates all 
renormalizations such as Pauli principle induced pair-pair re- 
pulsion, pairing symmetry and density related effects. Note, 
in contrast to Ref. |^, in the present work we are not restricted 
to the bosonic limit, nor is it essential to consider a periodic 
lattice. Thus, much of this language is also relevant to the 
moderately strong coupling (but still fermionic) regime, and 
can even be applied to jellium. 

The goal of this subsection is to establish a natural frame- 
work for relating M* to Tc. The parameters which enter into 
M* via Eq. ( |l2[ ) vary according to the length scales in the var- 
ious physical models. In the case of jellium, M* depends in 
an important way on the ratio k^/kp- For the case of s-wave 
pairing on a lattice, M* depends on the inverse lattice constant 
n/a and density n. Finally, for the case of d-wave pairing. 



there is an additional length scale introduced as a result of the 
finite spatial extent of the pair. This enters as if there were an 
equivalent reduction in fco / in the analogous j ellium model . 
The following factors act to increase M* or, alternatively, to 
reduce the mobility of the pairs: the presence of a periodic lat- 
tice, a spatially extended pairing symmetry (such as d-wave) 
or, for jellium, small values of the ratio ko/kp < 0.4 (i.e., 
high density). 

In order to relate T^ to M*, we observe that in an ideal 
Bose-Einstein system Tc is inversely proportional to the mass. 
Here, this dependence is maintained, in a much more complex 
theory, as a consequence of Eqs. ( pa| ) and (p^). This is essen- 
tially an equation for the number of pairs (bosons), with renor- 
malized mass M*. Thus, as we increase g towards the bosonic 
regime, it is not surprising that Tc varies inversely with M*. 

This leads to our main observations, which apply to moder- 
ate and large g, although not necessarily in the strict bosonic 
regime, (i) For the general lattice case, we find that Tc van- 
ishes, either asymptotically or abruptly, as the coupling in- 
creases, in the same way that the inverse pair mass approaches 
zero.Lj (ii) For the case of jellium or low densities on a lattice, 
both Tc and M* remain finite and are inversely proportional. 
These observations are consistent with, but go beyond, the 
physical picture in Ref. ^ that Tc is expected to be propor- 
tional to the pair hopping integral t'. It should be stressed that 
in the very weak coupling limit the pair size or correlation 
length is large. In this case, the motion of the pairs becomes 
highly collective, so that the effective pair mass is very small. 

In the presence of a lattice, the dependence on band filling 
n is also important for M*, and thereby, for Tc- We find that 
the bosonic regime is not accessed for large n > nc ~ 0.53. 
There are two reasons why superconductivity abruptly disap- 
pears within the fermionic regime. This occurs primarily (in 
the language of Ref. ||) as a consequence of large pair-pair re- 
pulsion, relevant for high electronic densities, which leads to 
large M* . In addition, there are effects associated with the 
particle-hole symmetry at half filling.!^ Precisely at half fill- 
ing (i.e., the "filling factor" / = 1/2, or 2/ = n = 1), for 
the band structure we consider, there is complete particle-hole 
symmetry and /i is pinned at Ep- Similarly, in the vicinity of 
n = 1, the chemical potential remains near Ep for very large 
coupling constants g. 

By contrast, in the small density (lattice) limit for the s- 
wave case, (n 0.1), pair-pair repulsion is relatively unim- 
portant in M* and there is no particle-hole symmetry. In this 
way the bosonic regime is readily accessed. Moreover, in this 
limit we see a precise scaling of Tc with 1 /g in the same way 
as predicted by Ref. || (via the parameter t' = —2t^/U). Thus 
in this low density limit superconductivity disappears asymp- 
totically, rather than abruptly. 

The effects of pairing symmetry should also be stressed. 
Because of the spatial extent of the d-wave function, the pair 
mobility is strongly suppressed, and, thus, M* is relatively 
larger than for the s-wave case. This lower mobility of d- 
wave pairs leads to the important result that superconductivity 
is always abruptly (rather than asymptotically) destroyed with 
sufficiently large coupling. Near half filling we find /i remains 
large when Tc vanishes, at large g. As the density n is reduced. 
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FIG. 1. Tc and ra/M* as a function of p/gc in the 3d jellium 
model with ko/kp = 4 (main figure) and ko/kp = 1/3 (inset), 
corresponding to short range (or low density), and long range inter- 
actions (or high density), respectively. 

away from half filling, /i decreases somewhat. It is important 
to note that the system remains in the fermionic regime (with 
positive II) for all densities down to n « 0.09. 

In all cases discussed thus far, Tc exhibits a nonmonotonic 
dependence on the coupling constant. It grows exponentially 
at small g and shuts off either asymptotically or abruptly at 
higher g. One can view this effect as deriving from a com- 
petition between pairing energy scales and effective mass or 
mobility energy scales. This competition is not entirely dis- 
similar to that found in more conventional Eliashberg theory 
where the fermionic renormalized mass and the attractive in- 
teraction compete in such a way as to lead to a saturation in Tc 
at large coupling. However, in the present context, for inter- 
mediate and strong coupling, we are far from the Fermi liquid 
regime and the effective mass of the quasibound or bound pair 
is a more appropriate variable. 

With this background, it should not be surprising that non- 
monotonic behavior will arise, even in situations as sim- 
ple as in jellium models. Indeed, in this case we find 
that for sufficiently long range interactions or high densities 
(small ko/kp) superconductivity disappears abruptly before 
the bosonic regime can be reachedEj Even for the case of 
short range interactions {kg/kp — 4), there is a depression in 
Tc caused by an increase in the pair mass, while still in the 
fermionic regime. 

In Fig. |l] we plot the calculated Tc for the case of an 
isotropic, 3d jellium model with s-wave pairing, along with 
the inverse pair mass ra/M*. This figure is presented pri- 
marily as a base line with which to compare subsequent plots. 
The parameter k^/kp — 4, is reasonably large so that the 
high g asymptote is found to reach the ideal Bose-Einstein 
limit (Tc = 0.218£;f) with M* = 2m. The apnioach to 
the high g asymptote is from below, as is expected.Q This is 
a result of the decreasing Pauli principle repulsion associated 
with increasing g, and concomitant reduction in pair size. The 
nonmonotonic behavior at intermediate g/gc « lean be asso- 
ciated with structure in the effective pair mass, apd has been 
discussed previously from a different perspective.!!^ 

In the inset are plotted analogous curves for the case of long 



range interactions or high densities (ko/kp = 1/3). This fig- 
ure illustrates how superconductivity vanishes abruptly before 
the bosonicjegime is reached, as a consequence of a diverging 
pair massE3 

B. Effects of dimensionality 

In this subsection we illustrate the effects of anisotropy or 
dimensionality on Tc (and on Apg and /i ) within the context 
of a jellium dispersion.L^ A particularly important check on 
our theoretical interpolation scheme is to ascertain that Tc is 
zero in the strict 2d limit and that ii varies continuously from 
Ep in weak coupling to the large negative values character- 
istic of the strong coupling bosonic limit. The present cal- 
culational scheme should be compared with that of Yamada 
and co-workersO who included "mode coupling" or feedback 
contributions to Tc, but only at the level of the lowest order 
"box" diagram discussed in Ref. |l2[ These authors were un- 
able to find a smooth interpolation between weak and-strong 
coupling, but did successfully repair the problemsEjIla asso- 
ciated with the NSR scheme, which led to negative /i even in 
arbitrarily weak coupling. 

Figures ^(a) andg(b) show the effect on Tc and on Apj, and 
/i, respectively, of introducing a layering or anisotropy into 
jellium with s-wave pairing. The various curves correspond 
to different values of the anisotropy ratio TO^/my. It can be 
seen from these two figures that Tc approaches zero as the 
dimensionality approaches 2. At the same time the chemical 
potential /i interpolates smoothly from the Fermi energy at 
weak coupling towards zero at around g/go = 1.5 to large 
negative values (not shown) at even larger g. The vanishing of 
the superconducting transition in strictly 2d was discussed in 
detail in Sec. 0. 

It should be noted that quasi-two dimensionality will be an 
important feature as we begin to incorporate the complexity of 
rf-wave pairing. The essential physics introduced by decreas- 
ing the dimensionality is the reduction in energy scales for 
Tc. The chemical potential and pseudogap amplitudje-are rel- 
atively unaffected by dimensional crossover effects.EJ While 
Tc rapidly falls off when anisotropy is first introduced into 
a 3d system (such as is plotted in Fig. the approach to 
the strict 2d limit is logarithmic and therefore slow, as can be 
seen explicitly in Fig. ||(c). Thus, in this regime, to get fur- 
ther significant reductions in Tc associated with a dimension- 
ality reduction requires extremely large changes in the mass 
anisotropy. 



C. Effects of a periodic lattice 

The first appUcations of a BCS Bose-Einstein crossover the- 
ory to a periodic lattice were presented in Ref. ^ The present 
approach represents an extension of the NSR theory in two 
important ways: we introduce mode coupling or full self- 
energy effects which are parametrized by Apg, and which en- 
ter via Eq. (^. Moreover, the number equation [see Eq. (p3). 
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FIG. 2. Dimensionality crossover in a quasi-2djellium model, (a) 
Tc as a function of g is seen to vanish for all g as m±/m\\ —* oo, 
while (b) fi and Apg change little. A continuous variation of 
Tc versus m±/m^^ at g/go = 4 is shown in the main por- 
tion (semi-log plot) and the inset (log-log plot) of (c). Here 
ko/kp = 4, go = — 47r/mfco. 

which is a rewriting of Eq. (Q) in terms of Apg] is evalu- 
ated by including self-energy effects to all orders. This is in 
contrast to the approximate number equation used in Ref. ^ 
which includes only the first order correction. In this way we 
are able to capture the effects which were qualitatively treated 
by these authors and which are associated with the lattice. 

Figure |](a) plots the behavior of Tc (solid line) in 
an isotropic three-dimensional lattice (with s-wave pairing, 
(fik = 1) at a low density n = 0.1. The effects of higher 




-g/6t {g/gj 

FIG. 3. (a) Tc and (b) m/M* (solid lines) vs g at low filling 
(n = 0.1) on a 3d lattice, and Tc at larger filling in the inset of (a). 
A fit to the functional form t' = —2t^/g is plotted (dashed lines) 
in (a) and (b) with adjusted proportionality constants. For compari- 
son, m/M* vs g/gc for 3d jellium (Fig.ll|) is replotted (dotted line) 
in (b). From bottom to top, the inset of (a) shows Tc for densities 
n = 0.2, 0.5, 0.7, 0.85, and 0.9. The inset of (b) shows at = 
as a function of n. 

electronic filling are shown in the inset. The low n behav- 
ior in the main portion of the figure can be compared with 
the jellium calculations of Fig. |l]. For small n, Tc decreases 
asymptotically to zero at high g. For larger n, Tc vanishes 
abruptly before the bosonic regime (/i < 0) is reached [See 
inset of Fig. ||(b)]. These various effects reflect the analogous 
reduction in the effective pair mobility, parametrized by the 
inverse pair mass m/M* . To see the coiTelation with m/M* 
in the low density limit, we plot this quantity in Fig. 
for the lattice as well as jellium case (where for the latter, 
m/M* 1/2 at large g). Here the coupling constants are in- 
dicated in terms of g/gc for jellium and —g/6t for the lattice. 
The inflection points at —g/Qt « 2 in both Tc and m/M* 
curves correspond to = 0, which marks the onset of the 
bosonic regime. 

Also plotted in both Fig. p](a) and||(b) (dashed lines) is the 
effective hopping t' — —2t /g for n = 0.1, rescaled such that 
it coincides with Tc and m/M* , respectively, at high coupling 
(—g/Qt = 30). This figure illustrates clearly the effect first 
noted by Nozieres and Schmitt-Rink that in the entire bosonic 
regime, Tc varies with high precision as t' or equivalently as 
m/M*. 

Finally, in the inset of Fig. |^(b), we demonstrate the lim- 
iting value of n, above which the bosonic limit can not be 
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FIG. 4. Lattice effects on Tc (main figure) and jj, and Apg (inset) 
as a function of g for n=0.7 (dashed lines), 0.85 (dotted lines), and 
0.9 (solid lines) in quasi-2d for (a) s-wave and (b) d-wave pairing 
symmetries. Here tx/fy = 0.01. In (b), Tc vanishes at a much 
smaller g than does its s-wave counterpart. 

accessed. What is plotted here is the value of Tc at which /i 
is zero as a function of density n. This figure indicates that 
the bosonic regime can not be reached for n > ric ~ 0.53. 
At densities higher than this, the pair-pair repulsion increases 
M* sufficiently, so that vanishes abruptly, while is still 
positive. 

D. Effects of rf-wave symmetry 

We now introduce the effects of a rf-wave pairing interac- 
tion. For the purposes of comparison we begin by illustrating 
Tc for the case of s-wave pairing on an anisotropic lattice, 
shown in Fig. Q(a), fer three different values (0.7, 0.85, and 
0.9) of the density n.E3 The inset indicates the behavior of the 
pseudogap magnitude and the chemical potential. The plots 
of Apg for the three different n are essentially unresolvable in 
the figure. Note, from the inset, that within small numerical 
errors Tc and vanish simultaneously. A comparison of the 
magnitude of Tc (in the main figure) with the 3d counterpart 
shown in the inset of Fig. ||(alillustrates how Tc is suppressed 
by quasi-two dimensionality.Ej 

In Fig. 0(b), similar plots are presented for the d-wave case. 
Here we use the same values of the filling factor as in Fig. ia). 



to which Fig. ^b) should be compared. The essential differ- 
ence between the two figures is the large g behavior. Lattice 
effects produce the expected cutoff for s-wave pairing. In the 
d-wave situation this cutoff is at even smaller g, and moreover, 
corresponds to /i « Ep- Calculations similar to those shown 
in the inset of Fig. ^(b) indicate that superconductivity disap- 
pears while fi remains positive for all n ahpve the extreme low 
density limit (i.e., for n > Uc ~ 0.09)Ej This behavior is in 
contrast to that of the s-wave case where Uc « 0.53. 

In the rf-wave case, the pair size cannot be made arbitrarily 
small, no matter how strong the interaction. As a result of 
the extended size of the pairs, residual repulsive interactions 
play a more important role. In this way, the pair mobility is 
reduced and the pair mass increased. Thus, as a consequence 
of the finite pair size, in the d-wave case the system essentially 
never reaches the superconducting bosonic regime. 

E. Phase diagrams 

In this section we introduce an additional energy scale T*, 
and in this way, arrive at plots of characteristic "phase dia- 
grams" for the crossover problem. Our focus is on the pseudo- 
gap onset, so that attention is restricted to relatively small and 
intermediate coupling constants g; consequently, the bosonic 
regime is not addressed. Here, our calculations of T* are 
based on the solution of Eq. (pb|), along with Eq. (pc|), under 
the assumption that Apg = 0. This approxiruation for T* is 
consistent with more detailed numerical workEJ in which this 
temperature is associated with the onset of a pair resonance in 
the T-matrix. 

In Figs. |(a)^(c) our results are consolidated into phase 
diagrams for the different physical situations. The case of 
3d jellium, with s-wave pairing [Fig. ^(a)] is presented pri- 
marily as a point of comparison. Figure ^(b) corresponds to 
quasi-2d jellium {m^/m\\ = 10''), with s-wave pairing and 
Fig. ^(c) to the case pf rf-wave pairing in a quasi-2d lattice 
case (i±/i|| = 10^'').EZI The insets indicate the behavior of /i 
and Apg. Comparing T* with Tc represents a convenient way 
of determining the onset of the pseudogap state. (For definite- 
ness, we define the onset to correspond to T* — l.lTc). It is 
clear from the first two figures that this occurs for 3d jelliiim 
at g/ gc ~ 0.9, and for the quasi-2d case at g/go ~ 0.4.E3 
This observation reinforces the notion that pseudogap effects 
are easier to come by in lower-dimensional systems. Similar 
behavior is seen in the quasi-2d lattice situation for the d-wave 
case, although the energy scales on the horizontal and vertical 
axes reflect the parameter <|| (rather than gc and Ep). 

IV. IMPLICATIONS FOR THE CUPRATES 

There has been much concern in the literature about 
whether generalized BCS Bose-Einstein crossover theories 
are relevant to the copper oxide superconductors. Is the cou- 
pling g sufficiently "large" in some sense to warrant this form 
of departure from conventional BCS theory? In more concrete 
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FIG. 5. Phase diagrams of (a) 3d jellium (ko/kp ~ 4), quasi-2d 
jellium (ko/kp = 4,m±/m\^ = 10'*), (c) a d-wave symmetry on a 
quasi-2d lattice (t±/t^\ = 10"*). Here we take n — 0.9. The same 
energy units are (a), (b) Ep and (c) 4t y . 

terms, one may ask if the calculated energy scales for Apg, /i, 
T*, and Tc are consistent with experiment? Are there other 
effects which are more important than is the role of small ^? 
Perhaps among the most intriguing questions raised is how 
does one incorporate hole concentration (denoted by x) de- 
pendences into this picture? 

An early motivation for adopting these crossover ap- 
proaches was the observed short coherence length ^, which 
was suggestive of some form of "real space pairing." -It is 
also clear that these systems are doped Mott insulatorsEJ so 
that the metal insulator transition at 1/2 filling (x = 0) should 
be integrated.into any theoretical approach. This transition 
is generallyCJO parametrized through an "order parameter" 
such as the plasma frequency uip which must necessarily van- 




FIG. 6. Doping dependence of Tc, T* , and Apg [= 2 Apg, the mag- 
nitude of the pseudogap at (tt, 0)]. Here tx/ty = 0.01, to = 0.5 eV, 
and g is assumed constant, while the symmetry is d-wave. 

ish as a; ^ 0. Finally, it should bejipted that there are no 
dramatic effects on ^ with variable xJt3 

It is clear that, in any attempt to understand pseudogap phe- 
nomena in the cuprates, both the small size of ^ and that of ujp 
shouldJ^e addressed on an equal footing. Early work by our 
groupcj investigated the effects of small ujp on the crossover 
problem, at the level of the Nozieres Schmitt-Rink approx- 
imation, for charged fermions. Coulomb interactions were 
treated in the RPA in parallel with the RPA-like ladder dia- 
grams of the particle-particle attraction. These early calcu- 
lations established that deviations from the BCS limit were 
more pronounced, the smaller the plasma frequency. Thus, 
for the same value of g, proximity to the insulating state is 
correlated with a tendency towards "bosonic" superconduct- 
ing transitions. Additional effects (in the same direction) re- 
sult from the likely decrease in dimensionality as the insulator 
is approached. Both of these effects, thus, suggest an amplifi- 
cation of pseudogap phenomena as x 0. 

In this section we address two issues which have been 
raised as relevant for the cuprates: we examine the size of 
the various energy scales for the case of d-wave superconduc- 
tors (all of which depend on the hopping matrix element 
and we discuss some aspects of the hole concentration depen- 
dence of the phase diagram, beyond the very general qualita- 
tive issues which have been noted above. It should be stressed 
that the fundamental basis of all crossover theories is mecha- 
nism independent. No information is used or assumed about 
the details of the pairing mechanism beyond the existence of 
a pairing coupling constant g. In almost any microscopically 
based pairing scenario, g is likely to contain some degree of 
X dependence. However, since there is no consensus on the 
pairing mechanism, in the present paper it is inappropriate to 
obscure our general results by making any detailed assump- 
tions about the nature of g{x). 

Here we focus exclusively on the x dependence of the un- 
derlying metal-insulator transition. We take g as doping in- 
dependent (which is not unreasonable in the absence of any 
more detailed information) and incorporate the Mott transi- 
tion at half filling, by introducing an a;-dependence into the 
in-plane hopping matrix elements i|| of our calculations. In 
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this way, we can explore the question of the size of the vari- 
ous energy scales and capture some degree of hole concentra- 
tion dependence, albeit not the entire effect. Our renormalized 
band structure is based on the limit of extremely strong on-site 
Coulomb repulsion, as seems appropriate for representing the 
Mott transition. It follows from very early work on the Hub- 
bard modeO that the hopping matrix element is renormalized 
as <|| (x) « to{l — n) ~ Iqx, where to(~ 0.5eV) is the ma- 
trix element in the absence of Coulomb effects. Equivalently, 
the effective particle mass varies as 1 /x. This change of en- 
ergy scale is consistent with the requirement that the plasma 
frequency vanish at x = O.M 

In Fig. ^ we replot the d-wave phase diagram of Fig. ^ 
for the case of fixed ~g/4to = 0.045 (and t^/t\\ = 0.01), 
which is chosen to fit the measured size, of the pseudogap at 
Tc for extremely underdoped cuprates.Ea Shown in the figure 
are T*, Apg, and Tc- Agreement with experiment may or may 
not be fortuitous since the coupling constant was assumed to 
be independent of x. Nevertheless, the energy scales, appear 
to be consistent with those measured experipientallycjLlJ and 
the X dependent trends are not inconsistent.^ 

It should be stressed that the results shown in the figure are 
robust consequences of our crossover theory. As a result of 
d-wave symmetry, Tc vanishes at moderately strong coupling. 
Moreover, this maximal coupling is a fairly universal number 
(i.e., independent of n) for a given over the physical range 
of hole concentrations (x < 0.3). [See, e.g.. Fig. H(b)]. Once 
a linear x dependence is enforced in t\\, near half filling, Tc 
decreases naturally as x decreases, and vanishes for extremely 
low doping concentration. This feature is insensitive to the 
detailed parametrizations of the model. 



V. CONCLUSIONS 

In this paper we have applied a previously discussedl^ 
BCS Bose-Einstein crossover theory to complex situations 
which are more physically relevant than are our earlier stud- 
ies of 3d s-wave jellium. In this way, we have determined the 
effects of quasi-two dimensionality, of periodic discrete lat- 
tices, and of a d-wave pairing interaction. This crossover the- 
ory yields results which appear consistent with known phys- 
ical constraints and plausibility arguments. Thus, in partic- 
ular, our strict 2d calculations yield a sensible interpolation 
scheme (for ^) with Tc strictly zero. The effects of the lattice 
are consistent with earlier Monte Carlo and other approaches, 
yielding a vanishing strong coupling limit for Tc associated 
with an increase (with g) of the effective mass M* of the 
fermion pairs. Finally our d-wave studies reveal that, for this 
symmetry, the superconducting bosonic regime is essentially 
never reached. Tc is suppressed to zero at moderate coupling 
constants, presumably because of the lowered pair mobility 
due to the constraints imposed by d-wave symmetry: the pair 
size cannot be reduced beyond the scale of a lattice spacing. 
These features should be appended to other observations in 
the literatures which note that in the strong coupling limit a 
d-wave superconductor will not exhibit gap nodes. Indeed, 



it is sometimes argued that this provides a "proof that the 
cuprates (which exhibit explicit d-wave symmetry) cannot be 
in the bosonic regime. Our results appear to make this case 
even more strongly, since we find that Tc will be zero when- 
ever a d-wave system is in the preformed pair limit. 

Our paper includes a brief discussion of the relevance to the 
copper oxide superconductors, wherein we impose the sim- 
plest possible ingredients of a Mott transition to arrive at some 
indications of hole concentration dependence and characteris- 
tic energy scale parameters, such as Tc, T*, and the pseudogap 
amplitude Apg. The numbers which emerge seem to be rea- 
sonably consistent with experiment, although we have made 
no assumptions about the origin or hole concentration depen- 
dence of the pairing interaction. In this way, one may argue 
that these crossover scenarios provide useful insights into the 
pseudogap state of the cuprates. 
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